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Abstract. We propose hypergeometric constructions of simultaneous approximations 
to polylogarithms. These approximations suit for computing the values of polyloga- 
rithms and satisfy 4-term Apery-like (polynomial) recursions. 



The series for the logarithm function 



is so simple and nice that mathematicians immediately generalize it by introducing 
the polylogarithms 



and considering then multiple, g-basic, p-adic and any other possible generalizations, 
just to have a serious mathematical research (i.e. to have some fun). It is so that 
nobody could now overview the whole amount of the results around all these gener- 
alizations of the logarithm, since the literature on the subject increases to infinity as 
a geometric progression (almost hypergeometrically) . 

It is not surprising that the transcendence number theory also dreams of getting 
new and new results for the values of the polylogarithms, especially after Lindemann's 
proof of the transcendence of log x for any algebraic x different from and 1 . The 
main problems (or, if you like, intrigues) are therefore extensions of the result to the 
values of (1), where Lindemann's argument based on proving the transcendence for 
the inverse, exponential, function does not work in an obvious manner any more. 
And we now have, so far irrationality and linear independence results for the poly- 
logarithm values at nonzero rational points close to zero, thanks to contributions of 
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Maier [Ma] (1927), Galoclikin [Ga] (1974), Nikisliin [Ni] (1978), Chudnovsky [Ch] 
(1979), Hata [Hal], [Ha2] (1990 and 1993), Rhin and Viola [RV] (2004). The very last 
piece of news is the irrationality of Li2(l/(?) for q integer, q ^ —5 or q ^ Q, obtained by 
the powerful arithmetic method in [RV], which improves the range of [Ha2] by adding 
q = 6 (the work [RV] also includes quantitative improvements of the irrationality in 
other cases, but we do not touch this subject in this short note). Another direction 
of arithmetic investigations are the values of (1) at 2; = 1 (or 2; = —1), so called zeta 
values. This goes back to Euler's time, who has contributed by the formula 

Li,,(l) = C(2/c) = i^il^^^M!!^ for k = l,2,..., (2) 

where i?2fc G Q are the Bernoulli numbers, thus Lindemann's proof of the transcen- 
dence of TT results in the transcendence of the numbers (2). Apery [Ap] (1978) has 
shown that C(3) is irrational, and since that time, thanks to Ball and Rivoal [BR], we 
dispose of only partial irrationality information for other values of Li2fc-i-i (1) = C(2A;-|-1) 
if = 2,3,4, ... . 

All known achievements in this subject are closely related to hypergeometric series 
and also multiple and complex integrals originated from the series. This is a general 
concept of the hypergeometric method developed for arithmetic study of the values 
of the polylogarithms; we refer the reader to a brief exposition of this concept in [Ne]. 

Here we would like to present some new ingredients of the hypergeometric method. 
We cannot achieve some new number-theoretic results by these means, and for the 
moment this note may be viewed as a methodological contribution. Nevertheless, 
approximations to the values of the polylogarithms that we derive here are quite 
reasonable from the computational point of view, and, in this sense, we continue our 
previous work on deducing curious Apery-like recurrences. 

We hope that the reader is somehow familiar with our work on the hypergeometric 
method in arithmetic study of zeta values (at least with the preprint [Zul]). 

1. Simultaneous approximations to the logarithm and dilogarithm 
For each n = 0, 1, . . . , consider the rational function 

((t-l)(t-2)---(t-n))2 



Rnit) = 



nl-t{t + l)---{t + n) 



Since degree of the numerator is greater than degree of the denominator, we will 
have a polynomial part while decomposing into partial fractions. The arithmetic 
properties of this decomposition are given in the following statement; denotes the 
least common multiple of the numbers 1, 2, . . . , n. 

Lemma 1. We have 

- n\-t{t + l)---{t + n) -^TTk^ 

fe=0 

where numbers are all integers and ■ B{t) is an integer-valued polynomial of 
degree n — 1. 
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Proof. Write this decomposition as follows: 



n-1 



t{t+l)---{t + j-l) 



t + k 



(the empty product for j = is 1). The coefficients Ak are easily determined by the 
standard procedure: 



Ak^Rnit){t + k 



t=-k 



/c = 0,f,...,n, (3) 



while the remaining group of unknown coefficients requires some work. Denote 

Ak 



Fiit) = it + l).J2 



k=0 



t + k'' 



I = 0,1, . . . ,n. 



Then 



and 



dFijt) 
dt 



k=0 
k^l 



k-l 



Therefore, 



dt 



t=-i k=0 
ki^l 



k-l 



What will happen if we do the same with the polynomial tail? Define 



n-l 



Gi{t) = {t + l)-Y,Bi 



t{t+l)---{t + j -1) 



l = 0,l,...,n-l. 



Then 



dGiit) 



dt 



t=-l 

and, since for a polynomial P{t) 

d 



i=o 



tit+l)---{t + j-l) 



dt»— ^ + 



t=-l 



dt 



we finally obtain 

dGijt) 
dt 



t=-i 



{Pirn + 1)) 



t=-l 



B 



I = 0,1, . . . ,n — 1, 



3=0 



(4) 
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and hence 



Bi 



dGi{t) 



dt 



(5) 



Furthermore, we have 



Fi{t)+Gi{t)^Rn{t){t + l) 



it-l)---it-n) 



n!(t + 



t{t+l)---{t + ny 



Z = 0, 1, . . . , n — 1, 



hence 



Dn--{Fi{t)+Gi{t)) 



e Z, / = 0, l,...,n-l, 



t=-l 



where Dn denotes the least common multiple of the numbers 1, 2, . . . , n. These inclu- 
sions and inclusions 



Dr.- 



dFijt) 
dt 



e Z, / = 0,l,...,n-l, 



t=-l 



which follow from formulae (3) and (4), together with the induction on I on the basis 
of (5), show that 

Dn-BieZ, / = 0,l,...,n-l, 



and the proof follows. 

Since for an integer- valued polynomial P{t) of degree at most n its derivative 
multiplied by Dn is again an integer-valued polynomial, we also have 

Lemma 2. The following decomposition is valid: 



dRnjt) 

dt 



fe=0 



A, 



{t + ky 



+ B{t), 



where numbers Ak are all integers and ■ B{t) is an integer-valued polynomial of 
degree n — 2. 

Let 2; be a rational number with < \z\ < 1. We are now interested in the following 
two hypergeometric-type series: 



rn = rn{z) = ^ z"" Rn{t)\^^^, r^ = rn{z) = - 



' dt 
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Lemma 3. We have 

rn{z) = anUi{z) - bn, rn{z) = anU2{z) -bn, 

Z^an e Z, iziZ2yDnbn G Z, {ziZ2TDlbn G Z, 

where z\ and zi are the denominators of the numbers 1/ z and zj (1 — z^, respectively. 
Proof. Let us write the polynomials B{t) and B{t) in the form 

where _ 

Dn-BjeZ, Dl-B^eZ, j = 0,...,n-l (7) 

(this is guaranteed by the theorem of choosing a basis in the Z-space of integer- valued 
polynomials). Then 

fc=0 1^=1 j=0 v=l 



k \ I n—1 



Z^+' 



E'— — E 

fc=0 fc=0 1=1 j=0 

In the same vein, 

fc=0 fe=0 i=l j=0 

Using (7) and integrality of for all A; = 0, 1, ... , n, we arrive at the desired claim. 



fc=0 fc=0 1=1 j=0 ^ ^ 

ne vein, 

" " ^-{k-l) '^-2 _ / „ \ : 

; — n ; — n ; — i — n \ / 



Remark. As follows from (3) and the above proof, we have the following explicit 
formula: 

El n\ f n + k\ ( 1 



k=0 



kj \ k J \ z 



As we see from Lemma 3, the sequences r^iz) and r^(-z), n = 0,1,..., realize 
simultaneous rational approximations to Lii(2) and Li2(-2), although approximation 
'tails' are not simply polynomials in 1/z, but sums of two polynomials from Q[l/-2;] 
and Q[z/{z — 1)]. Although we worked in the region \z\ < 1, the result of Lemma 3 
remains valid in the closed disc l-^l ^ 1 except the point 2; = 1 by analytic continuation. 
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Running the Gosper-Zeilberger algorithm of creative telescoping [PWZ, Chapter 6] 
with the input Rn{t)z*, one can find the difference operator of order 3, which anni- 
hilates linear forms (6) and their coeflBcients an,bn,bn- In order not to frighten the 
reader, we indicate only the characteristic polynomial of this operator 

z{z - 1)X^ - {3z^ - 20z + 16)A2 + z{3z + 8)A - z^ 

(containing all details on the asymptotic behaviour of the approximants, due to 
Poincare's theorem), and the partial case 2; = — 1 of the corresponding recurrence. 

Theorem 1. The coefficients an, bn and bn of the simultaneous approximations 

rn = rn{-l) = anLii(-l) -bn = -an log2 - bn, 

TT^ ~ n = 0, 1,..., 

Tn = fn{-'^) = a-n Ll2(-1) - &n = - 

as well as the approximations themselves satisfy the recurrence relation 

2(59n - 24) (n + ifa^+i ~ (2301n^ + 136571^ - 376?i - 240)a„ 

- (295n^ - 120n2 - 60n + 35)a^_i - (59n + 35)(n - 1)^0^-2 = 0, 

n = 2, 3, . . . , 

of order 3, and the necessary initial data is as follows: 

ao = 1, ai = 5, a2 — 55, 

n L 7 ^ 305 r n r . r 181 
bo = 0, bi = --, 62 = — oo = 0, 61 = -4, 62 = — ^• 



In addition, 



lim |r^|^/''= lim jF^I^/'^ = |Ai 2 H 0.15960248 ... , 



lim lonl^/'' = lim Ibnl^/"" = lim Ibnl^/"" = A3 = 19.62866250 ... , 

n—^oo n— >oo n^oo 

where Ai^2 = —0.06433125 . . . ± iO. 14606314 . . . and A3 are zeros of the characteristic 
polynomial 2X^ — 39A^ — 5A — 1. 

One can run the algorithm of creative telescoping with the input Rn{t) (that is, 
2; = 1, a non-sense!) to obtain a much simpler recurrence 

(n + l)^an+i -I- (lln^ -|- lln -|- 3)an — n^an-i = 0, n = 1, 2, . . . , 

of order 2, which may be recognized in view of Apery's proof of the irrationality 
of C(2) (see [Ap] and [Po]). In fact, we have the identity 
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thanks to Thomae's transformation of 3^2(1 )-hypergeometric series (see [Ba, Sec- 
tion 3.2]), and the latter sum gives (up to the sign factor) the denominators of Apery's 
approximations to C(2)- In order to give the necessary sense to the substitution z = 1, 
we should introduce the following complex Barnes integral: 

TT^ / Rn{t)i-. -] z'dt = r n{z)-rn{z) log z, (8) 

27rzic-ioo \sm7rtj 

where C is an arbitrary constant in the interval < C < n+1. The integral converges 
in the whole disc \z\ ^ 1, hence we should have the limit in the right-hand side of (8), 
i.e. (in notations of the proof of Lemma 3) the limit 

exists and is equal to a certain rational constant depending on n. On the other hand, 
the complex integral in (8) admits a real double-integral representation thanks to [Ne, 
Theorem 2]. Taking m = 3, r = 2 and ai = a2 = = n + 1, b2 = = 2n + 2 in this 
Nesterenko's theorem, we obtain 

rn{z)-Uz)logZ = jj (i_^ + ,^^)n+l ^xdy. 

[0,1]2 

Substituting z = 1 and multiplying by (—1)" reduces the latter integral to Beukers' 
famous double integral [Be] for Apery's approximations to C(2)- 

2. Simultaneous approximations to C(2) and ({3) 

Our first natural generalization of the construction in the previous section is based 
on the rational function 

P _ {{t-l){t-2)...{t-n)r 
nP-t{t + l)---{t + n) ■ 

Then 'reasonable' approximations to the first three polylogarithms are given by the 
series 



'^n{z) = ^z''Rnit)\^^^ = anUi{z) - h„, 



z 



.w = ^E 



i/=i 

oo 

n /-(■„ I r I 

= anl-'hiz) - hr. 

oo ,2 

r 1 I / 1 

= an Lis (2;) - br. 



i/=i 



1^ ,d^Rn{t) 



2^ dt2 



where 

Z^ttn e Z, {ziZ2)'^Dnbn G Z, {ziZ2)'^DnD2 eZ, {zrZ2TDnDlj}n&'^, (9) 
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and zi and Z2 are the denominators of the numbers 1/z and z/{l — z), respectively. 
The reason of having the multiples in (9) is the higher degree 2n — 1 of the 
polynomial in the decomposition of Rn{t) into partial fractions, and it is required to 
derivate it for getting the representation of rvi and rn- The explicit formula for the 
coefficient a„ is as follows: 



i~rt \ / \ 



k I \ z 



k=0 

This time we are interested in the particular 'non-sense' case 2; = 1 of the con- 
struction. Without tiring reader's eyes by writing two complex integrals converging 
in the disc |z| ^ 1 and guaranteeing the existence of limits for corresponding series 
expansions, we just present the final result for the approximation sequences 

rn = rn{l) = anCC^) -K, rn = rn{l) = anCi^) -K, n = 0, 1, 



n 



Theorem 2. The above sequences rn, Vn as well as the coefficients an, bn and h 
satisfy the recurrence relation 

2(946n^ - 731n + 153) (2n + 1) (n + ifan+i 

- 2(104060?!^ + 127710n^ + 12788?i'* - 34525?i^ - 848271^ + 3298n + 1071)an 
+ 2(3784?!^ - 1032?!^ - 1925?i^ + 853n'^ + 328n - 184)nan-i 

- (946n2 + 1161n + 368)n(n - 1)^0^-2 = 0, n = 2, 3, . . . , 
of order 3, and the necessary initial data is as follows: 

ao = 1, ai = 7, a2 = 163, 
~ ~ 23 ~ 2145 r . ~ 17 ~ 3135 

60 = 0, 61 = y , 62 = bo = 0, 61 = y , 62 = y^. 

In addition, 

limsup Irnl^/"" = limsup IfnT^" = |Ai,2| = 0.067442248 . . . , 

n — >oo n — >oo 

lim |an|^/" = lim = lim = A3 = 54.96369509 ... , 

n—K30 n— >oo n— >oo 

where Ai,2 = 0.018152450 . . .±zO. 064953409 . . . and A3 are zeros of the characteristic 
polynomial 4:X^ - 220X'^ + 8A - 1. 

Based on the recurrence, we have observed experimentally and we are able to show 
that the correct form of the inclusions (9) in this special case 2 = 1 is 

»-Mr|:MK:)Cr)'=|:(:)\"r)("r)-, 

DnD2nbn & Z, Olbn G Z. 

Remark. Normalizing the approximations of Theorem 3 by multiplying them by the 
factor (^^), we arrive at the recurrence previously obtained in [Zu2, Theorem 4] by 
means of a certain implicit construction. Our new explicit consideration gives an 
answer to arithmetic observations posed in [Zu2]. 
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3. Well-poised approximations 

The arithmetic study of zeta values was strongly influenced by well-poised hyper- 
geometric series. They are a 'heart' of the proof in [BR] and of many other similar 
results, and we cannot avoid considering a well-poised generalization of the construc- 
tion in Section 1. 

Take 

nm ( ^^n+l^f , n\ ((t - 1) • • • (t - n) • (t + n+ 1) • • • (t + 2n))^ 

= {-irRn{-t-n), 
which has now a quite complicated partial-fraction decomposition: 

where (repeating arguments of the proof of Lemma 1) 2Ak, 2D„^'^, 2D'^A'^ are inte- 
gers for all A; = 0, 1, . . . , n, while 2Df^B{t) is an integer-valued polynomial. Gathering 
this arithmetic knowledge and proceeding as in the proof of Lemma 3, we deduce that 
the series 



are certain linear forms involving certain polylogarithms (up to Li4(2;)). Not so 
exciting, but we would like to deal with the construction at the only one point, 
z — —1. The well-poised thread (take — t — n in place of t in (10)) results in equal- 
ities (-1)^=^^ = -(-1)— and A'^ = -{-1)^-''A';_,,, A; = 0, 1, . . . , n, 

and they are the circumstance, which makes r„(— 1) and rn{—l) linear forms in 
Li2(— 1) = — 7r^/12, 1 and 2Li3(— 1) = — 3C(3)/2, 1, respectively, with the same 
leading coefficient. We write this final production as follows: 

rn = rn{-l) = - bn, Tn = rn{-l) = an^^Y^ -K, 77, = 0, 1, . . . , 

where 

2DnaneZ, 2''DlbneZ, 2''D%eZ, n = 0,l,...; (11) 

the nth powers of 2 appear since the two is the denominator of z/{z — l) when z = —1. 
Applying the algorithm of creative telescoping we arrive at the following result. 

Theorem 3. The sequences Tn, Tn and the coefficients an, bn and bn satisfy the 
recurrence relation 

(1457n^ - 1363n + 348) (n + ifan+i 

- (148614n^ + 158202n^ - 9295n^ - 6189471^ - llllln^ + 8932n + 2784)an 
+ (97619n^ - 91321n^ - 9443n^ + 3534371^ - 5440/1^ - 5678n + 1768)a„_i 

- 3(1457n^ + 1551n + 442)(3n - 2)(3n - 4)(n - 1)^0^-2 = 0, n = 2, 3, . . . , 
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of order 3, and the initial values are as follows: 



ao = 1, ai = 8, a2 = 264, 

13 , 1737 ~ ~ 29 ~ 7617 

oo = 0, 6i = y , 62 = — 60 = 0, 61 = y , 62 = y^. 

In addition, 



limsup Irnl^/'' = limsup |fn|^/" = |Ai,2| = 0.51616460 ... , 

n— >oo n— >oo 

lim |on|^/"= lim l^n^^" = lim = A3 = 101.34149804 ... , 

n— >oo n^oo n— >oo 

where Ai^2 = 0.32925097 . . . ± *0. 39751691 . . . and A3 are 2;eros 0/ the characteristic 
polynomial A^ - 102A2 + 67A - 27. 

On the basis of the recurrence relation we find much better inclusions than (11), 
namely 

i-^a^eZ, n = 0,l,..., (12) 

where is the following product over primes: 



= n 

2/3^{n/p}<l 

{ ■ } denotes the fractional part of a number. Inclusions (12) are quite expected in view 
of 'denominator conjectures' around linear forms in zeta values (see [KR, Section 17.1] 
about the difficulties in proving the correct arithmetic in similar cases). But why do 
we get the cancellation of 2"'? This might be also caused by the well-poised origin of 
the series used by us. At least the integrality of a„ is an immediate consequence of 
the following explicit formulae: 




where equality (13) follows from [KR, Proposition 5] and equality (14) is communi- 
cated to us by G. Almkvist. 
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4. Final remarks 

As already promised, no new irrationality and linear independence results were 
presented. We just have tried to give some sense to certain hypergeometric-type 
series that are expressed in terms of polylogarithms and are divergent when one 
formally plugs z with 1^1 = 1. Transforming a non-terminating single hypergeometric 
series into a multiple one (some kind of 'identites non-terminees gigantesques', cf. 
[KR, Section 17.5]) often meets convergence troubles for the latter series, i.e. it is 
just a formal transformation, which we could never use in a rigorous proof. As an 
option to proceed in such troubling cases, we see dealing with transformations for 
complex Barnes (multiple) integrals and further decompositions of the integrals into 
sums involving (multiple) zeta values. This does not look an easy program, but we 
do not believe that deducing new results for zeta values and polylogarithms might be 
simple. 



References 

[Ap] R. Apery, Irrationalite de C(2) et C(3), Asterisque 61 (1979), 11-13. 

[Ba] W. N. Bailey, Generalized hypergeometric series, Cambridge Math. Tracts, vol. 32, Cambridge 
Univ. Press, Cambridge, 1935; 2nd reprinted edition, Stechert-Hafner, New York, 1964. 

[BR] K. Ball and T. Rivoal, Irrationalite d'une infinite de valeurs de la fonction zeta aux entiers 
impairs, Invent. Math. 146 (2001), no. 1, 193-207. 

[Be] F. Beukers, A note on the irrationality of Q(2) and C(3), Bull. London Math. Soc. 11 (1979), 
no. 3, 268-272. 

[Ch] G. V. Chudnovsky, Pade approximations to the generalized hypergeometric functions. I, J. 

Math. Pures Appl. (9) 58 (1979), no. 4, 445-476. 
[Ga] A. I. Galochkin, Lower estimates of linear forms in values of certain G-functions, Mat. Za- 

metki [Math. Notes] 18 (1975), no. 4, 541-552. 
[Hal] M. Hata, On the linear independence of the values of poly logarithmic functions, J. Math. 

Pures Appl. (9) 69 (1990), no. 2, 133-173. 
[Ha2] M. Hata, Rational approximations to the dilogarithm. Trans. Amer. Math. Soc. 336 (1993), 

no. 1, 363-387. 

[KR] C. Krattenthaler and T. Rivoal, Hypergeometric et fonction zeta de Riemann, Preprint (Feb- 



ruary 2004), E-print fnath. NT/0311 114| y3. 

[Ma] W. Maier, Potenzreihen irrationalen Grenzwertes, J. Reine Angew. Math. 156 (1927), 93— 
148. 

[Ne] Yu. V. Nesterenko, Integral identities and constructions of approximations to zeta-values, J. 

Theorie Nombres Bordeaux 15 (2003), no. 2, 535-550. 
[Ni] E. M. Nikishin, On irrationality of values of functions F(x,s), Mat. Sb. [Russian Acad. Sci. 

Sb. Math.] 109 (1979), no. 3, 410-417. 
[PWZ] M. Petkovsek, H. S. Will and D. Zeilberger, A = B, A. K. Peters, Ltd., Wellesley, MA, 1996. 
[Po] A. van der Poorten, A proof that Euler missed... Apery's proof of the irrationality of C(3), 

An informal report. Math. Intelligencer 1 (1978/79), no. 4, 195-203. 
[RV] G. Rhin and C. Viola, The permutation group method for the dilogarithm, Prepublication 

no. 2004-13, Departement de Math., Metz, 2004. 
[Zul] W. Zudilin, An elementary proof of Apery's theorem. Preprint (February 2002), E-print 



nath. NT/0202159 



Zu2] W. Zudilin, A third-order Ap ery-like recursion for C(5), Mat. Zametki [Math. Notes] 72 



(2002), no. 5, 733-737; E-print [math. NT/0206178 



Department of Mechanics and Mathematics, Moscow Lomonosov State University, 
VoROBiovY Gory, GSP-2, 119992 Moscow, RUSSIA 
E-mail address: wadimSips.ras.ru 



